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ABSTRACT 

We investigate the color-magnitude (CM) relation of galaxies in the distant X-ray selected cluster 
' RDCS 1252.9-2927 at z=1.24 using images obtained with the Advanced Camera for Surveys (ACS) on the 

I Hubble Space Telescope in the F775W and F850LP bandpasses. We select galaxies based on morphological 

. classifications extending about 3.5 mag down the galaxy luminosity function, augmented by spectroscopic mem- 

' bership information. At the core of the cluster is an extensive early-type galaxy population surrounding a central 

' pair of galaxies that show signs of dynamical interaction. The early-type population defines a tight sequence in 

^ ■ the CM diagram, with an intrinsic scatter in observed (/775-Z85o) of 0.029 ±0.007 mag based on 52 galaxies, or 

' 0.024 ±0.008 mag for ~30 ellipticals. Simulations using the latest stellar population models indicate an age 

scatter for the ellipticals of about 34%, with a mean age rj, > 2.6 Gyr (corresponding to zl Si 2.7), and the last 
star formation occurring at Zend ^1-5. Transforming to rest-frame (U-B), we conclude that the slope and scatter 
I I in the CM relation for morphologically selected early-type galaxies show little or no evidence for evolution out to 

^ . z « L2. Thus, elliptical galaxies were akeady well established in X-ray luminous clusters when the universe was 

\^ ' a third of its present age. 

■ Subject headings: galaxies: clusters: individual (RDCS 1252.9-2927) — galaxies: elliptical and lenticular, cD 
I — galaxies: fundamental parameters — cosmology: observations 

o : 

, 1. INTRODUCTION galaxy clusters with spectroscopic confirmation. This Letter 

■ T, . , 1 . II- 1 1 ui 11 u I, J presents the first results from our ACS cluster survey, focusing 
Present-day cluster ellipticals are a remarkably well-behaved ^ , , . , , , . ^ , , ■' , ° 

^ , f u- » -.u . . 1 J u • 1 u on the color-magnitude (CM) relation of the early-type galax- 

,C class of obiects, with structural and chemical properties obey- . . „,^„„,^?^ , , , , 

r^ . • 1 1 1- w D . .u- ij . 1 les in RDCS 1252. We adopt the best-fit WMAP cosmology: 

ing simple power-law scaling relations. But this could not al- rm A^^^ /t> . i iAna\ ■ • 

I • u u »u ■ u- u- 1 • \\ru■^ » (/I, S 2,,,, 0/,) = (0.71 , 0.27, 0.73) (Bennett et al. 2003), givmg a 

A I ways have been the case in a hierarchical universe. While most ^ \ ^^r, 

^ I f J 1 U.J. J scale of 8.4 kpc per arcsec at z= 1.237. 

^ galaxy formation models can be tuned to reproduce these re- ^ ^ 

lations at z=0, a more stringent test lies in reproducing their 

. evolution with redshift. To this end, it is important to study 2. OBSERVATIONS AND IMAGE REDUCTIONS 

' rich clusters out to the highest redshifts, when fractional age DT>>r-c nci u a ■ ,u tmc^r jcocnrou a 

■ j.r-j, .1, 1 ■ . 1 . T RDCS 1252 was observed in the F775W and F850LP band- 
. differences among the galaxies were proportionately greater In ., ^ . , . , s . , , . , 

, jjcij.i A A A- passes (hereafter i775 and Z850, respectively) with the ACS Wide 

^r' • recent years, deep wide-field optical surveys and deep serendip- i^. , , ^ ^ ■ ■ . , 

.J , J VC. C r u Field Camera as part of the guaranteed time observation pro- 

?H itous X-ray surveys have uncovered significant numbers of rich , , „i,„„, , . ^„„^ , , , ^^^^ , rr^, 

rri ' 1 1 . . J u f. J u A / ■ u gram (proposal 9290) during 2002 May and 2002 June. The 

'■^ , galaxy clusters to redshift unity and beyond (see reviews by °, ^t- t- , . ^ ^ . . , - , ^ 

S , onm r. onm\ t-u . j- . . j . obscrvations Were done in a 2 x2 mosaic pattern. With 3 and 5 

Postman 2002; Rosati 2003). These most distant, and most ,. ^. ... , ., 

c, 11 u J . . .u u orbits of integration in i77g and Z850, respectively, at each of the 
massive, of known gravitationally bound structures can then be ^ . . ° ^, , ' , , \ 
. J- J • J . 1 .u u . . J u- u 1 f 11 four pointings. There was nearly 1 of overlap between point- 
studied in detail through targeted, high-resolution, follow-up . ^, ,° ^ . -'. ■j- ,^,^ ,. 

... .cji- .- ings; thus, the core of cluster was imaged for a total of 12 orbits 

optical and near-infrared observations. .^^^ 2q Qj^j^g jjj 

We have undertaken a survey of rich galaxy clusters in the ''JJ^ , 

juvc. no^ ^ 1 o ■ .uAj A c The data wcrc proccsscd With the Apsis pipeline described 
redshift range 0.8 < z < 13 using the Advanced Camera for , , , , /^^^^x . , , 
o /A^c T- J . 1 onm\ .u rr ; ; ; c i by Blakcslcc ct al. (2003), With somc rcccnt Updates. In par- 
Surveys (ACS; Ford et al. 2002) on the Hubble Space Tele- / , , ■ r , , ■ , n , c, 
/rrc"TA -ru ■ f .u- • . . ui- u ticular, wc used a version of the drizzle software (Fruchter & 
scope (HST). The aim of this survey is to establish new con- , ,. , , ,t . . . , . 

'. , .u 1 . f u J .u 1 f Hook 2002) supplied by R. Hook that implements the Lanc- 

straints on the cluster formation epoch and the evolution of -„ . , • , w i i • ^ ■ n • i 

, ^ , . c . 1 . 1- J ni^^o 1 '-ic'-i n zos3 interpolation kcmcl (a damped sinc function). This ker- 

early-type galaxies. The first cluster observed, RDCS 1252.9- , , ^ , . ■ /^.r^r-N i i 

2927 (hereafter RDCS 1252) at z= 1 .237 (Rosati 2003; Rosati P'""'^'!'^^^ ^. ^hai-per point spread function (PSF) and greatly 

. 1 -irvmN J- J . f.u r./-.cAT^i^ /-^i . reduccs the noisc coiTclation of ad] acent pixcls and the resulting 

et al. 2003), was discovered as part of the ROSAT Deep Cluster . „ ... . .. . , 

e /r. . 1 inno\ j • .u u- u f j u vc. moirc patterns. Apsis also now rcmovcs discontinuities in the 

Survey (Rosati et al. 1998) and is among the highest-redshift , , , , ■ , ■ 

residual bias level at the amplifier boundaries, producing a more 
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Fig. 1 . — Color composite of the core region of RDCS 1252, constructed from our ACSA¥FC F775W and F850LP images, shown in the observed orientation. 
The displayed field size is roughly 1' across, or less than 4% of the full mosaic. 
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uniform background. An earlier processing of these images has 
been used by Bouwens et al. (2003) for a study of the faint (775 
dropout population at z ~ 6. We calibrate our photometry to 
the AB system using photometric zero points of 25.640 (1775) 
and 24.843 (zsso)- These are uncertain at the ^^0.02 mag level, 
which has no effect on our conclusions. We adopt a Galac- 
tic reddening for this field of E{i—z) = 0.041 mag based on the 
Schlegel et al. (1998) dust maps. 

3. OBJECT SELECTION AND PHOTOMETRY 

Figure 1 shows the central -^1' region of a color composite 
made from our reduced i^^5 and zgso images. A red galaxy pop- 
ulation is clearly visible. The central pair of galaxies are sep- 
arated by l."8 (15kpc) and are each of magnitude Z850 ~ 21. 
We used SExtractor (Bertin & Arnouts 1996) in "dual-image 
mode" with low threshold and deblending settings to find ob- 
jects in the reduced images and perform the initial photometry. 
SExtractor "MAG_AUTO" values were used for the total mag- 
nitudes. The (/775-Z85o) color effectively separates out evolved 
galaxies at z > 1, and the cluster is obvious as a central con- 
centration of galaxies with 0.80 < (j775-Z85o) < 1-05. Figure lb 
[removed from the ApJL version in order to meet the page limit] 
shows histograms of isophotal color within 3 different radii 
of the cluster center (defined midway between the two central 
galaxies) for galaxies with total Z850 = 20-25 mag. 

We selected an initial sample of 312 nonstellar objects 
with Z850 < 24.8, in the broad isophotal color range 0.5 < 
('775-285o) < 1-2, and inside a radius of 1.'92. Our goal is 
to study the early-type galaxy population in RDCS 1252, for 
which we have limited spectroscopic data, and these cuts are de- 
signed to select the vast majority of our target sample while rea- 
sonably limiting foreground/background contamination. The 
color selection is roughly 7 times broader than the full- width of 
the red sequence we find below. The radial cutoff corresponds 
to about 1 .0 Mpc for both our adopted WMAP cosmology and 
an Einstein-deSitter cosmology with h=Q.5. 
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"Fig. lb." — Histograms of (1775-Z850) color for objects having 20 < zgso < 
25 and within the specified radii of the center of RDCS 1252. The cluster galax- 
ies are near (;775-Z85o) - 0.95, and make up an increasing fraction of the objects 
within progressively smaller radii. [This figure was removed from the published 
version because of space limitations.] 



Our final colors are measured within galaxy effective radii Rg 
to avoid biasing the CM slope due to color gradients. We follow 
the basic approach outlined by van Dokkum et al. (1998, 2000). 
We derive the Re values using the program "galfit" (Peng et al. 
2002) by fitting each galaxy to a Sersic model (convolved with 
the PSF), but constraining the n parameter such that 1 < n < 4. 
Bright neighboring galaxies were fitted simultaneously. We 
note that subtraction of the model for the two central galax- 
ies reveals evidence for interaction in the form of an S-shaped 
residual. 

Next, we deconvolve (775 and Z850 postage stamp images of 
each galaxy using the CLEAN algorithm (Hogbom 1974) in or- 
der to remove the differential blurring effects of the PSF, which 
is ~ 10% broader in the Z85o band. To reduce noise, the CLEAN 
maps are smoothed with a Gaussian of FWHM= 1.5 pix before 
adding the residual images to the maps to ensure flux conser- 
vation. We then measure the iJ^5 and zsso magnitudes of each 
galaxy within a circular aperture of radius R^, typically about 6 
pix, or 0"3. We did not allow the radius to drop below 3 pix. 
Photometric errors were determined empirically from the point- 
ing overlap regions. We reprocessed the 4 pointings separately 
and measured the color differences within Re for 202 pairs of 
measurements for 74 different early-type galaxies (classifica- 
tions described below) in the overlap regions. We then median- 
filtered to obtain the errors as a smooth function of magnitude. 
The error thus determined for single-pointing (ins—z&so) mea- 
surements at Z850 = 23 was 0.025 mag, rising to ~ 0.05 mag at 
Z850 = 24.5. 

Each galaxy in our initial sample was examined and mor- 
phologically classified, following a procedure similar to that of 
Fabricant et al. (2000). This was done independently of the pro- 
file fitting, but the types show a good correlation with Sersic n 
index. Here, we simply distinguish between E, SO, and later 
types, where the intermediate SO class indicates the apparent 
presence of a disk without spiral or other structure. Full de- 
tails on the classifications for a much larger sample, including 
the spatial distribution of the various types, will be presented 
by Postman et al. (2003, in preparation). About 180 galaxies 
in this field have measured redshifts, obtained with VLT/FORS 
(Rosati et al. 2003, in preparation), with 31 being cluster mem- 
bers. All galaxies in our initial sample classified as early-type, 
and not known to be interlopers from their spectra, are included 
in our CM analysis in the following section. Of the 3 1 known 
members, 22 are classified as early-type, and we include all of 
these in our analysis even though one (an SO) happened to lie 
beyond our 1 Mpc cutoff. 



4. THE COLOR-MAGNITUDE RELATION AT z= 1 .24 

We fitted the early-type galaxy CM relations using simple 
linear least squares; other methods gave very similar results. 
We estimate the scatter from both the standard rms and the bi- 
weight scale estimator (Beers et al. 1990). No rejection was 
done in fitting subsamples composed of known members, the 
faintest of which has Z85() = 23.48, or ^0.5L^. We also per- 
formed fits to samples with unconfirmed members, allowing us 
to go 1 mag further down the galaxy luminosity function. How- 
ever, here we iterate to reject the 3-(t outliers, as these are likely 
to be interlopers: none of the 22 confirmed early-type members 
is more than 2.3-a discordant. After the iterative rejection pro- 
cess, we find concordant scatters for those samples, and the rms 
and biweight estimator are the same to within ±0.001 mag. 
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Figure 2 presents the CM relation for the RDCS 1252 galax- 
ies. The fit to the full sample of early-type galaxies with 
Z850 < 24.5 gives 

(i-z) = (0.958 ±0.006)- (0.025 ±0.006)(z85o- 23). (1) 
Other results are listed in Table 1 . The mean locations of the 
CM relations for the elliptical and SO subsamples agree to well 
within the errors, while the slopes are consistent at the l.S-a 
level. Eight known late-type members from Rosati et al. (2003, 
in preparation) for which we have photometry are bluer than 
the early-type galaxies by 0.25 mag, with a scatter of 0.14 mag 
about this offset, indicating young stellar populations. We find 
an intrinsic scatter (Tint = 0.023 ±0.007 mag for the 15 con- 
firmed elliptical members. For a Umit of zsso < 24.5, we de- 
rive (Tint = 0.026 mag for the clipped sample of 31 ellipticals 
and CTint = 0.029 for the 52 E+SO galaxies. At this limit, the 
observational errors become dominant and classification is dif- 
ficult, which could bias our tJint estimates. For zs5o < 24.0, still 
3 mag down the luminosity function to about 0.3 L*, we find 
fint = 0.024 for 25 ellipticals (with no outliers). 
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Fig. 2. — Color-magnitude diagram for early-type galaxies inside 1^9 with 
O775-Z850) > 0.5, excluding spectroscopicaUy known interlopers. Circles and 
squares represent ellipticals and SOs, respectively; solid symbols are used in 
the CM relation fits, while open symbols (all of which lack spectroscopic in- 
formation) are rejected as outliers or as below the faint cutoff (indicated by 
the dotted line). Two representative fits are shown: the fit to the 15 ellipti- 
cal members (solid line) and to the 52 early-type red-sequence galaxies. The 
approximate luminosity conversion for RDCS 1252 is shown at top, assuming 
the WMAP cosmology and —1.4 mag of luminosity evolution as described in 
the text, such that Mg = -21.7 (AB). The relation for the Coma cluster, trans- 
formed to these bandpasses at z = 1.24 (no evolution correction), is indicated 
by the dot-dashed line. 

We estimate « -21.7 AB for RDCS 1252, based on lo- 
cal surveys (Norberg et al. 2002), WMAP cosmology, and -1.4 
mag luminosity evolution (Postman et al. 2001; van Dokkum 
& Stanford 2003). The correction from Z850 to rest-frame B is 
-1-0.3 mag for an early-type spectrum. Thus the brightest clus- 
ter member at zsso = 21.0 corresponds to ~4.8L^. The clus- 
ter red-sequence is about 0.1 mag bluer than predicted for non- 
evolving elliptical templates (Coleman et al. 1980; Schneider 
et al. 1983). For comparison, a Bruzual & Chariot (2003; here- 
after BC03) solar metalhcity model reddens by 0.11 mag and 
fades by 1 .5 mag in aging from ~ 3 to ~ 11 .5 Gyr. The present- 
day relation for the Coma cluster (transformed according to the 
relations given below) is also shown in Figure 2. 

The tight CM relation allows us to constrain the scatter in 
early-type galaxy ages, subject to model uncertainties. The evo- 
lution in (j775-Z85o) at this redshift is complicated at young ages 



because Z850 straddles the 4000 A break and has its blue end 
near the Balmer jump at 3700 A. The Balmer jump reaches a 
maximum at 0.5-1 Gyr, when the relative contribution from A 
stars is greatest, resulting in a red color at these ages. It is inter- 
esting that several of the SOs lie above the CM relation, possibly 
indicating ages < 1 Gyr. As the A star contribution lessens, the 
color quickly evolves towards the blue, reaching a local mini- 
mum near ~ 1.5 Gyr, before commencing a roughly monotonic 
reddening with age. 

We have simulated the (z'tts-Zsso) colors of galaxies formed 
under two simple star formation models, similar to those used 
by van Dokkum et al. (1998). In Model 1, the galaxies form 
in single bursts randomly distributed over the interval ito,hnd), 
where fo is set to the recombination epoch and fend < f^ = 5.1 
Gyr, the age of the universe at z = 1.24. In Model 2, galaxies 
form stars at constant rates between randomly selected times 
iti,t2), where fo < ^1 < ^2 < ?end- We vary fend and at each step 
calculate colors and luminosities for 10,000 "galaxies" by inter- 
polation and integration of the BC03 solar metallicity models. 
Assuming amt = 0.024 mag for the ellipticals. Model 1 imples 
a minimum age fz-fend =1-6 Gyr, i.e., all galaxies finish form- 
ing at redshifts z > Zend = 1-9; the mean luminosity-weighted 
age is T£ = 3.3 Gyr (corresponding to zl = 3.6) with a scatter of 
30%. Model 2 gives f^-fend = 0.53 Gyr, Zend = 1 -4, and tl = 2.6 
Gyr (corresponding to zl = 2.7) with 38% scatter. Thus, al- 
though some galaxies in Model 2 have formed stars recently, 
the mean ages are still high. Both models give a mean color 
('775-Z85o) = 0.94, similar to that observed. 

For the SOs, we find Zend =1-5 for Model 1 and Zend = 1-3, 
indicating recent star formation, and age scatters of 44-47%. 
Finally, we note that the 1996 version of the BC models would 
have predicted higher formation epochs, e.g., Zend > 2.5 and 
Zend > 2.0 for the ellipticals in Models 1 and 2, respectively, and 
an age scatter of only 20%, although the predicted color is then 
redder by 0.1 mag. Overall, we conclude that the ellipticals are 
an evolved population, with a mean age tl ^ 2.6, a minimum 
age ~ 1 ±0.5 Gyr, and an age scatter of (34 ± 15)%, where the 
error reflects uncertainty in (Tint and scatter in the models. 



5. DISCUSSION 

To enable comparison with previous work, we convert ob- 
served (j775-Z85o) quantities to rest-frame (U-B)^ at z = 1 .24. 
The models (BC03; Kodama et al. 1998) and empirical tem- 
plates indicate A(U-B), = (1 .8 ± 0.4) x A(/775-Z85o), where the 
error bar reflects the scatter in the models at the relevant ages 
and adds about 20% uncertainty to our transformed slope and 
scatter. Figure 3 uses this conversion, and other transforma- 
tions from van Dokkum et al. (2000), to compare our results 
to some previous studies of the CM relation in intermediate- 
redshift clusters with HST, as well as the results of van Dokkum 
et al. 2001 on RX J0848H-4453, flie only cluster of comparable 
redshift to have its CM relation measured. 

Linear fits to the data shown in Figure 3 yield slopes of 
-0.014 ±0.010 and 0.003 ±0.008 for the evolution in the abso- 
lute slope and the scatter, respectively. Thus, the scatter is con- 
stant, and there is at best marginal evidence for slope evolution, 
which indicates that the slope is due to a variation in metallicity, 
not age. Previous studies of cluster samples out to z ^ 1 have 
come to similar conclusions (e.g., Stanford et al. 1998; Kodama 
et al. 1998; van Dokkum et al. 2000). Van Dokkum et al. (2001) 
concluded that the slope at z = 1.27 was shaUower than in the 
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Coma cluster. However, as shown in the figure, our slope mea- 
surement is consistent within the errors with both Coma and 
RXJ0848H-4453. Further studies of a diversity of clusters at 
similar redshifts are needed to explore this issue. 

The lack of evolution in the CM relation scatter can be ex- 
plained by progenitor bias (e.g., van Dokkum & Franx 2001): 
galaxies selected as early-type at any epoch will have old stel- 
lar populations, while the later-type progenitors of the youngest 
ellipticals today will not be selected. The result is an underes- 
timate in the color scatter for the progenitors of modern ellipti- 
cals, and thus overestimated ages. An upper hmit on the scatter 
for elliptical progenitors may be estimated from a fit to all con- 
firmed RDCS 1252 members; the result is 3^ times larger than 
for the early-type galaxies. A detailed study of the morpholog- 
ical fractions in RDCS 1252 (Postman et al. 2003, in prepara- 
tion) should help illuminate the magnitude of this bias. We also 
note that the two central ellipticals themselves, based on their 
proximity and irregular isophotes, appear likely to undergo dis- 
sipationless merger to form a single ~9L| galaxy, similar to 
local cD galaxies. 

We conclude that massive, evolved early-type galaxies were 
already present in rich clusters at z = 1 .24. Our simple models 
imply mean luminosity-weighted ages of 2.6-3.3 Gyr, corre- 
sponding to formation at z = 2.7-3.6. However, the (ins-zsso) 
color is not ideal for the redshift of RDCS 1252, being better 
suited to measuring the 4000A break at z ?a 1.1. Combining 
our ACS data with deep, high-resolution near-lR imaging of 
this field (Lidman et al. 2003) will enable a more robust assess- 
ment of early-type galaxy ages. In addition, further studies of 
the CM relations and morphologies of galaxies in other z ^ 1 
clusters are needed to improve the constraints on the formation 
epoch of cluster galaxies and on the evolution of their stellar 
populations and structural properties. 
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Fig. 3. — Slope (top) and scatter (bottom) of the rest-frame (U—B) CM re- 
lation as a function of redshift. The filled circle is from the present work; open 
symbols show results from, in order of increasing redshift, Bower et al. (1991); 
van Dokkum et al. (1998); EUis et al. (1998); van Dokkum et al. (2000); and 
van Dokkum et al. (2001). The dotted lines indicate the average values. 
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Table 1 

RDCS 1252-2927 Color-Magnitude Relations 



Sample zf^^ Nt Nc Slope CTbwt" fbwt'' CTmt 







15 


15 


-0.022 ±0.011 


0.029 


0.029 ±0.005 


0.023 ±0.007 


E+SO". . 




22 


22 


-0.018 ±0.013 


0.038 


0.038 ±0.006 


0.033 ±0.007 


E" 


24.0 


25 


25 


-0.020 ±0.009 


0.033 


0.033 ±0.005 


0.024 ±0.008 


E+SO'! . 


24.0 


44 


41 


-0.025 ±0.008 


0.045 


0.038 ±0.004 


0.029 ±0.007 


E" 


24.5 


34 


31 


-0.019 ±0.007 


0.053 


0.036 ±0.005 


0.026 ±0.008 


E+SOl . 


24.5 


58 


52 


-0.025 ±0.006 


0.054 


0.039 ±0.004 


0.029 ±0.007 


SO". . . . 


24.5 


24 


21 


-0.042 ±0.013 


0.058 


0.039 ±0.006 


0.032 ±0.008 



"Biweight scatter based on raw number of galaxies iVr. 
''Biweight scatter based on iVc galaxies after 3-a clipping. 
■^Spectroscopically confirmed members of specified type only. 

"All red-sequence objects (known interlopers omitted) of specified type, Z850 < Zgso^ and 
within the area of analysis. 



